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Abstract 

Phase structure and phase transitions in dense QCD are studied using the 
Cornwall- Jackiw-Tomboulis (CJT) potential in the improved ladder approximation. 
The gap function, the condensation energy and the structure of Cooper pairs are in- 
vestigated at finite temperature and density. Due to strong coupling effects at low 
densities, the gap, critical temperature and their ratio deviate from the weak coupling 
values. It is shown that the internal structure of Cooper pairs are robust against the 
thermal effects, despite the fact that the pairs are strongly correlated near the critical 
temperature. Also, the effect of the strange quark mass, m s on the phase diagram is ex- 
amined using a simple kinematical criterion. We discuss the behavior of the unlocking 
line, on which the CFL turns into the 2SC, through the variation of m s . 
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§1. Introduction 



The color superconducting phases at high baryon density have attracted much interest 
in high density QCD. In particular, the two types of the pairing pattern, namely color-flavor 
locked (CFL) pairing and iso-scalar (2SC) pairing have been studied by many authors and 



In the chiral limit, weak coupling analyses of the gap equations reveal the CFL dominance 
over the 2SC at zero temperature and the coincidence of critical temperatures of transitions 
to the QGP phase from the 2SC state and from the CFL state. The former is attributed 
to a competition between the number of degrees of freedom participating in the pairing 
correlation and the sizes of the gaps. The latter is due to two facts: The rapid vanishing 
of the pairing in the color sextet channel as the critical temperature T c is approached, and 
the disappearance of the nonlinearity of the gap equation in the diquark condensate near T c . 
However, in the low density regime, the sizes of gaps are modified from the weak coupling 
values owing to various strong coupling effect sP Also, the pairing in the color symmetric 
channel may become relevant at low densities, although it is suppressed by one coupling 
constant relative to the anti-triplet channel in the weak coupling regime. 

We examine in this article the Schwinger-Dyson (SD) approactP in the improved ladder 
approximation to the pairing phenomena at finite temperature and density. This approach, 
the SD equation in the Landau gauge with an improved running coupling constant, can 
reproduce the reasonable vacuum physics related to the chiral symmetry breaking as well 
as the perturbative result for the gap in the weak coupling regime. Various strong coupling 
effects manifest themselves for realistic densities in this approachP More complete and 
detailed analysis was done in RefP 

§2. Gaps, correlation and thermodynamic grand potential 

The gap equation is nothing but the self-consistency condition for the proper self energy 
A±(po,p) in the Nambu-Gor'kov base \P = (g,i7 2 7°g*) 



Here, is the in-medium gluon propagator, S is the full quark propagator, and i" 1 ^ = 
7 At A a /2 is the quark-gluon vertex. The gap matrix is defined by the off-diagonal component 
of the self energy S. The 2SC ansatz for the gap matrix is 



rich physics in these phases has been revealed so farP 



<£(Po, p) = ig 2 






(2-2) 
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Z\ + (_) is the positive (negative) energy component of the gap matrix A(po, p) = (Z\_ + A + ) — 
(A + — Z\_)7°7 • p. The color-flavor locked gap matrix is 

A C ± F L ( P0 ,P) = \5«8)Al{p ,p) + (s*$ - A 8 ± ( Po ,p). (2-3) 

Index (a, b, ■ ■ ■ ) and • • • ) represents color and flavor, respectively. We use the quasi- 
static approximation of the hard dense loop propagator in the Landau gaugeP^ Following 
Ref.P we replace the coupling g 2 with the momentum dependent effective coupling g 2 (p,q) 
using the Higashijima-Miransky prescription. 

fip ' q) = 1 ^ H(PiJ+f c )/A 2 y ^ = max (^)' ( 2 - 4 ) 

where (3q = (lliV c — 2Nf)/3, and pi plays the role of an infrared regulator. We adopted A = 
400 MeV and p 2 c = 1.5 A 2 for our numerical calculations. The correlation function describes 
the off-diagonal long range order in the system: 

Ml) = tanh (j(q-ri* + Al(q)/2T) (2-5) 



W(q-v) 2 + Ai(q) 

A(q) is the gap for the quasi-quark having momentum q = \q\. 4>+{p) are defined by 
(2SC|a"(p)a^(— p)|2SC) = (r 2 A2)"j i 0+(g), where a is the annihilation operator for quark. 
We have omitted the helicity indices here. We can calculate the coherence length from 
the correlation function which characterizes the size of the quark Cooper pairP® Us- 
ing the method developed by Cornwall, Jackiw and TomboulisP the thermodynamic po- 
tential par unit volume relative to the normal Fermi gas up to 2-loop order is given by 
5f2((i,T) = {TrLog [SSq 1 ] — ^Tr^X 1 ] }. S is bare quark propagator in the Fermi sea. 
Note that this expression is only valid at the stationary point determined by solving the gap 
equation for S. 

§3. Numerical results and discussion 

We present the numerical solutions for the gap equations at fi = 1000 MeV corresponding 
to about 80 times the normal nuclear saturation density po = 0.17 fm~ 3 . We display the gap 
function evaluated on the quasi-quark mass shell A±(p) for the 2SC state Fig.|U(a), and the 
singlet and octet gap functions, Ai±{jp) and A 8 ±(p), for the CFL state (b). The shapes of 
these gap functions are all similar as functions of momentum, except for their magnitudes 
l^sG )! < 1^0°) I < l^i (p) I- m order to see the characteristics of the phase transitions in 
detail, we show the temperature dependence of the gaps at the Fermi surface Fig. [T] (c) and 
that of the CJT condensation energies for the 2SC and CFL states (d). The phase transitions 
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Fig. 1. (a) Gap functions at fj, = lOOOMeV for the 2SC state. The values (-\k\, A_(\k\)) are 
plotted in left half of the figure, (b) Gap functions for the CFL state, (c) The temperature 
dependence of the gaps at the Fermi surface, (d) The condensation energy as function of 
temperature. 



to the QGP phase both from the 2SC state and from the CFL phase are of 2nd order, as in 
the BCS theory with a contact interaction: The derivatives of the gaps seem to diverge as T c 
is approached in (c), and the thermodynamic bulk quantities in the superconducting phases 
connect continuously to those of the normal QGP phase (d). The size of Cooper pair £ c , the 
root mean square radius of the Cooper pair with wave function <f> + is displayed in Fig. [21(a). 
£ c is not affected significantly by changes in the temperature. This shows the fact that the 
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Fig. 2. (a) The size of 2SC Cooper pair as a function of temperature at fx = lOOOMeV. (b) The 
calculated QCD phase diagram from the improved Schwinger-Dyson approach. 



Cooper pairs are robust against the thermal fluctuation, although the correlation between 
diquarks diverge as the critical temperature is approached because A(q) — > as T — > T c — 0. 
This fact also suggests the existence of the incoherent Bose gas of tightly bound diquarks 
just above T c . 

Fig. |21 (b) shows our result for the QCD phase diagrams for m s = 0, 150, 200 and 
250 MeV. In the chiral limit, there is only one line (black bold dots) dividing the (/x, T) 
plane into the QGP and CFL phases by 2nd order transition. The critical temperature 
is increased as one goes towards low density. This is due to the strong coupling effects 
including participation of antiquarks to pairing correlation. It is surprising that the critical 
temperatures of transitions to the QGP phase from the 2SC and CFL states exactly coincides. 
The unlocking critical line on which the CFL state turns into 2SC state by 1st order transition 
appears for finite value of m s , and shifts from the triangles to the squares, and to crosses, with 
the increase of m s . These were estimated using a kinematical criterion for the unlocking.®*" 
The strength of the 1st order transition becomes weaker as the quark density is increased, 
and the critical end point at which the 1st order transition terminates is located at \i — oo. 
If the strange quark mass is lower than 200 MeV, then chirally broken vacuum phase might 
be continuously connected to the CFL phase without phase transition. 
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